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I. OBJECTIVE: 

The purpose of this publication is to illustrate data compiled 
by the University of California at Berkeley in a study of twelve 
spring and winter wheat segments in the Dakotas, Kansas, and Oklahoma. 

II. BACKGROUND: 

A. The UC3 Study 

The study concerns itself with variations in temporal -spectral 
characteristics for a given crop. The primary objective was to 
measure and quantitatively describe the year-to-year variation in 
crop temporal -spectral response patterns due to any year-to-year 
variation in the physical environment. 

In order to determine the yearly variation in crop temporal - 
spectral response, it was necessary to develop a measuring system 
which is acquisition independent. One way this can be accomplished 
is by estimating a continuous spectral response curve for a crop 
as a function of time. From such a curve, one can determine major 
spectral events (spectral biostages) for a crop (e.g. maximum 
amplitude of a vegetation indicator, date of first green canopy 
detection) and assess the year-to-year variation in such events. 

The function used to generate the continuous temporal plot is 

c _ ,31 J2 T 1 2 + B3 

r - T e 


Where F = the value of the green vegetation indicator (GRABS), T = day 
of the year, and Bl, B2, B3 are coefficients estimated by a least 
squares procedure. This particular function was chosen for its appro- 
priate shape and its adaptability by logarithmic transformation to 
estimation by least squares linear regression. However, estimation 
by least squares linear regression proved to be unsatisfactory 
because there is no obvious variance term to associate with a 
temporal -spectral value derived from such a segment-average curve. 


1. Hay, C.M., et. al . , Annual Progress Report , Space Sciences Laboratory, 

University of California (Berkeley), November 1979. 


In addition, the curve fit is extremely poor in some cases. For 
these reasons a nonlinear least squares program was developed to 
estimate the coefficients of the function. For each wheat field, 
the following variables were calculated from the equations generated 
by the nonliri®ar curve-fitting program. 

F 3 maximum GRABS value 

max 

T * day of maximum GRABS 

iTiaX 

* day of first detectability 

T 3 day of harvest 

a 

The corresponding average values for each segment-year were calculated 
by taking the unweighted arithmetic mean over all sampled fields 
within that segment-year. These values are represented by the curves 
in the accompanying graphs. 

B. The Data 

The data set was prepared from Landsat MSS data and the segments 
were selected on the basis of: 1) the number of years of Landsat 

and ground truth data available, 2) adequate within-year acquisition 
histories, 3) diversity in number of drought stress years, and 4) geo- 
graphical diversity. The MSS data was sun angle and haze corrected, 
and a green vegetation indicator was calculated for each pixel. Wheat 
fields in specific segment-year combinations were then randomly sampled. 

Ground truth data was available for all the segments in years 
'77 and '78, but of the segments plotted only 1851 had ground truth 
in '76; For those segment-years in which ground truth was not available, 
the wheat fields were identified by analyst interpretation. 

CALCULATION OF CURVES 

The curves illustrated in the graphs were derived using a pro- 
cedure common to the entire data set. The data consists of twenty 
three sets of the following variables. 

T1 3 date of frrst detectability 

FI 3 GRABS threshold for detectability (=2) 

T2 3 date of peak GRABS 
F2 3 peak GRABS 


T3 = date of harvest 
F3 = FI (=2) 

These variables are essentially x and y coordinates for three points 
on a curve. Using these values for T (day of year) and F (GRABS) 
along with the function 

F = t B 1 e 82 jZ + 83 

one can solve for a system of three equations and three unknowns 
using gaussian elimination. 

General Form: Ax = b 

ln(Tl)*Bl + (T1) 2 *B2 + B3 = ln(Fl) 
ln(T2)*Bl + (T2) 2 *B2 + 83 = ln(F2) 
ln(T3)*Bl + (T3) 2 *B2 + B3 = ln(F3) 

Example: For segment 1175, year '76 (see Table 1.2), the following 

were mean values over 15 fields: 

T1=40.87 T2=95. 55 T3=165.36 Fl=2.00 F2=18.00 F3=2.00 

The System to be solved is: 

3.71040*81 + 1670.35681*82 + 1.00000*B3 = 0.69315 

4.55965*B1 + 9129.80273*82 + 1.00000*B3 = 2.89037 

5.1 081 2*B1 + 27343. 92969*B2 + 1.00000*B3 = 0.69315 

The Solution Vector is: 

B1 = 4.95862 B2 = -0.00027 B3 = -17.25308 


An arbitrary number of values for F can now be generated for arbitrary 
values of T. 


F . = T.Bl e 82 T ,‘ . 83 ii= , !30 

Plotting F. versus T. yields the desired curve for this particular 
set of T1 , T2, T3, FI, F2, F3. This procedure can similarly be 
extended to plot the remaining data. It should be noted that the 
curves represent mean values for the fields within a given segment-year. 
Judging from the standard deviations (see Table 1.2) of the parameters, 
a profile curve for an individual field within a segment-year will likely 
vary considerably from the mean profile curve for. that segment-year. 


IV, CONCLUSIONS: 


UCB concluded that the curves show no significant correlation 
between the spectral variable F m _. w and any of the temporal variables. 
This result denies support to ideas asserting that the relative 
earl i ness or lateness of a crop's development affects the value of 
its spectral peak. However,* significant correlations were observed 
between T max and and between T max and T . Alternatively, the 
relative earliness (or lateness) of the GRABS peak was associated 
with the relative earl i ness (or lateness) of both first detectability 
and harvest. 

Another observation made by UCB contends that within any segment, 
the year-to-year variation of crop growth characteristics were generally 
more significant than was the field-to-field variation within the same 
year. In addition, the nature and magnitude of the field-to-field vari- 
ability was found to be constant in two thirds of the cases studied. 
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DESCRIPTION OF GRAPHS 

In referencing the various graphs, the title in the upper right hand 
comer of each graph will be used. 

I. SPRING WHEAT 

A. "Segment: 1602-1 636 

In these graphs (one for each segment), each curve represents crop 
growth characteristics for a particular year (for curve identification 
refer to the legend on each graph). The solution vector associated 
with each curve as well as the location of the segment are also given. 

B. "Spring Wheat Average" 

This is a curve representing the average of all the curves plotted 
in the seven graphs "Segment: 1602-1686", 

C. "All Segments**All Years" 

These are all the curves of the seven graphs "Segment: 1602-1686" 
plotted on one graph. 

D. "Year: 1975-1978" 

These curves are those in "Segment: 1602-1686" but are grouped Such 
as to compare segments by year. That is, holding the year constant 
across all the segments. 

E. "Segment Averages by Year" 

This graph represents an average of the 1975 curves across all the 
segments, an average of 1976 -curves across all the segments, etc. 

F. ' "Yearly Averages by Segment" 

The curves represent a computed average of the curves within each 
segment for all years in which data was available for the segment. 

II. WINTER WHEAT 

A complete parallel of the preceding graphs follows for winter wheat. 







■_k'. *'.f 
































WINTER WHEAT 




















